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ABSTRACT 

While it is generally accepted that Type la supernovae are the result of the explosion of a carbon- 
oxygen White Dwarf accreting mass in a binary system, the details of their genesis still elude us, 
and the nature of the binary companion is uncertain. Kasen (2010) points out that the presence of 
a non-degenerate companion in the progenitor system could leave an observable trace: a flux excess 
in the early rise portion of the lightcurve caused by the ejecta impact with the companion itself. 
This excess would be observable only under favorable viewing angles, and its intensity depends on 
the nature of the companion. We searched for the signature of a non-degenerate companion in three 
years of Supernova Legacy Survey data by generating synthetic lightcurves accounting for the effects 
of shocking and comparing true and synthetic time series with Kolmogorov-Smirnov tests. Our most 
constraining result comes from noting that the shocking effect is more prominent in rest-frame B than 
V band: we rule out a contribution from white dwarf-red giant binary systems to Type la supernova 
explosions greater than 10% at 2a, and than 20% at 3er level. 



1. INTRODUCTION 

Type la supernova (SN la) explosions are marvelous 
astrophysical tools, and they currently offer the most 
preci se way of constraining dark energy ()Sullivan et al.l 
1201 lh . Today, several thousand SNe la have been ob- 
served and theoretical models an d simulations are pro - 
gressing rapidly (see for example lAlmgren et alJl2010f ). 
and many aspects of SN la explosions can be repro- 
duced in detail. However, these cosmic explosions, stud- 
ied for decades, still are not fully understood. Partic- 
ularly, we lack a solid understanding of the progenitor 
systems. There is consensus that SNe la arise from the 
thermonuclear explosion of a carbon-oxygen (C-O) white 

fbianco@lcogt.net 

1 Department of Physics, University of California Santa Bar- 
bara, Mail Code 9530, Santa Barbara CA 93106-9530 

2 Las Cumbres Observatory Global Telescope Network, Inc. 
6740 Cortona Dr. Suite 102, Santa Barbara, CA 93117 

3 Department of Physics (Astrophysics), University of Oxford, 
DWB, Keble Road, Oxford, OX1 3RH, UK 

4 Center for Astrophysics and Space Astronomy, University of 
Colorado, 593 UCB, Boulder, CO, 80309-0593, USA 

5 Departments of Physics and Astronomy, University of Cali- 
fornia, Berkeley 

6 Department of Astronomy and Astrophysics, University of 
Toronto, 50 St. George Street, Toronto, ON, M5S 3H4, Canada 

7 LPNHE, Universite Pierre et Marie Curie Paris 6, Univer- 
site Paris Diderot Paris 7, CNRS-IN2P3, 4 Place Jussieu, 75252 
Paris Cedex 05, France 

8 Universite Paris 11, Orsay, F-91405, France 

9 CPPM, CNRS-IN2P3 and Universite Aix-Marseille II, Case 
907, 13288 Marseille Cedex 9, France 

10 INAF - Osservatorio Astronomico di Roma, via Frascati 33, 
00040 Monteporzio (RM), Italy 

11 Australian Astronomical Observatory, P.O. Box 296, Ep- 
ping, NSW 1710, Australia 

12 CEA, Centre de Saclay, Irfu/SPP, F-91191 Gif-sur-Yvette, 
France 

13 LBNL, 1 Cyclotron Rd, Berkeley, CA 9472 

14 Network Information Operations, DRDC Ottawa, 3701 
Carling Avenue, Ottawa, ON, K1A 0Z4, Canada 

15 Department of Physics and Astronomy, University of Vic- 
toria, PO Box 3055 STN CSC, Victoria BC, V8T 1M8, Canada 



dwarf (WD), but is the WD in a binary system with 
a main sequence (MS) or red giant (RG) star, accret- 
ing mass from the companion to approach the Chan- 
drasekhar limit (single degenerate -SD- scenario)? Or 
is the explosion caused by the merger of two WDs in 
a compact binary system (double degenerate -DD- sce- 
nario) ? Constraining the progenitor scenarios is key 
for learning the details of SN la explosion physics, and 
to improve our understanding of the effects of environ- 
ment on SN la explosions and thus of the systematics 
that still affect t he constraints on cosmology derived 
from SN surveys (iKessier et al.l 120091; IGuv et all l20ll 
IWood-Vasev et al.l l2007t for" a review of SN la cosmology 
sec IHowellll2010D 

No progenitor system of a SN la has yet been observed 
prior to explosion: these binary systems would be very 
faint and undetectable, at this time, in extra-galactic 
surveys. Population synthesis and environment stud- 
ies have not been able to firmly set constraints on the 
SN la progenitors. From the theoretical point of view, 
generating SNe la in the DD scenario presents difficul- 
ties. The mass transfer only successfully leads to a de- 
flagration if it occurs at a rate significantly slower than 
the Eddington limit, throu gh the formation of a thick 
disk (|Nomoto fe; IbenlH98 5L and even then fine tuning 
of various parameters might be needed (see iToutl 120071 
for a brief review, and the references therein). 

Some observational evidence might already disfavor SD 
progenitors. While the WD accretes mass from a com- 
panion in the SD scenario, the system should emit X- 
ray radiation for an extended period of time. Under 
the assumptions of continuous duty cycle, and that all 
SD progenitors would emit in the X-ray, the SN la rate 
is too high by over an order of magnitude compared to 
the num ber of X-ray sources observ ed in nearby elliptical 
galaxies ■: CiH'nnuv & Bogdan 20 1 0l), as well a s soft X-ray 
sources m our own galaxy (|Di Stefanoll2010al ). While this 
evidence can be used to set upper limits to SD progen- 
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itors, | Di Stcfano (2010b) points out that there are too 
few super soft X-ray sources, sources with energy typi- 
cally 10 to 100 eV and luminosity 10 -37 to 10~ 38 erg/s, 
to account for the Type la rate within the DD scenario 
as well, suggesting instead that super soft X-ray radia- 
tion may not always be emitted in nuclear-burning white 
dwarf systems and that it could be absorbed within the 
system itself. A few peculiar SNe la have been observed 
in the past few years to produce 56 Ni masses close to, 
or in excess of, the 1.4M Q theoretical limit for a WD 
mass: the Chandrasckhar limit. These apparent super- 
Chandrasckhar SNe la may origina te from the merger 
of two WDs, in the DP scenario dHowell et all 120061: 
ISilverman et al ll201lHYuan et alJ[20M ~ 
IKasenI ( 20101 ) - hereinafter K10 - explores the effect 



that a non-degenerate companion star would have on the 
observables of the explosion, and shows that in the SD 
scenario, the presence of a companion may manifest it- 
self in the early days after the SN explosion as a flux 
excess. As the cloud of SN ejecta expands, it collides 
with the companion. This impact shocks the expanding 
material creating a hole in the otherwise optically thick 
ejecta shell through which radiation can escape. An ex- 
cess flux is produced in the shocked gas, propagating in 
the direction of the observer, and it should be detectable 
when the geometry is favorable and the observer looks 
into the companion star. As the the equilibrium temper- 
ature of the shocked debris is inversely proportional to 
the distance from the WD center to the power of 3/4 (see 
Eq. 7 in K10), the flux excess is larger the larger the sep- 
aration from the companion, and assuming Roche lobe 
overflow, RG companions will leave the most prominent 
signature. The intensity of the feature is higher in bluer 
bands (see Section [3]). The effect only lasts a few days, 
completely vanishing by 10 days after explosion. 

With detailed early time lightcurves we may be able 
to identify the progenitor system of a particular SN ex- 
plosion, singling out events generated by red giant pro- 
genitors when seen from favorable angles. Unfortunately, 
early detailed SN lightcurves, with daily or so cadence, 
are still rare. New s urveys like the PanST ARRS Medium 
Deep Survey fPSl. iPastorello et al.ll2010h and the Palo- 
mar Transient Factory (PTF. lLaw et"aLl I2009D provide 
well sampled early SN la lightcurves, which might lead 
to the identification of progenitors in individual cases, as 
might early UV follow up. 

With the large collection of lightcurves provided 
by sur veys such as the Supernova Legacy Survey 
(SNLS, lAstier et all I2006D and the Sloan Digital Sky 
Survey CSDSS. lAbazaiian et al.1 120091) the companion 
scenarios can be constrained in a statistical fash- 
ion. The SDSS collaboration recently searched for ev- 
idence of an early flux excess due to shocking in high 
signal-to- noise ratio (S/N) spectroscopically confirmed 
SN la lightcurves from the SDSS-II survey. Finding no 
evidence of shocking-related excess in a subset of 108 
confir med SNe la with w ell observed early time behav- 
ior, [Haydcn e t al.llMlOal conclude that RG's cannot be 
the main channel for SN la explosions. Here we use con- 
firmed SNe la from the first three years of SNLS data 
to set an upper limit to the contribution of RGs to 
SN la progenitors. 

After describing the dataset used here, consisting of 
spectroscopically confirmed SNe la from the first three 



years of SNLS, and the processes used to standardize the 
data and generate composite lightcurves (Section [2]), in 
Section[3]we summarize the results presented in K10 and 
show our rendering of these models. We then describe 
the statistical tests that allowed us to derive constraints 
to the contribution from RG progenitors to SN la ex- 
plosions (Section H]). We also extended our analysis 
beyond the spectroscopically confirmed sample to in- 
clude photometrically selected SN la lightcurves, show- 
ing that lightcurves affected by shock were not rejected 
as SN la spectroscopic follow up candidates in SNLS be- 
cause of a selection bias, and that our conclusions extend 
to the photometrically selected SNe la (Section [SJ . Fi- 
nally, in Section [7] we summarize our conclusions and 
outline future work. 

2. SNLS DATA, THIRD YEAR 
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The d ataset used 
tail in iConlev et al.l (1201 ID . 
IGonza lcz-Gai tan et al.l (|2011l ) and 
We use data from the first 3 years of the SNLS. 
The SNLS is a rolling survey that gathered photo- 
metric data at the Canada-France-Hawaii- Telescope 
(CFHT). Two independent photometric pipelines, 
based in France and Canada, are used for the SNL S 
data reduction (|Bazin et al.l 120111 : iPerrett et al~ll2010( ). 
Here we use the photometry output of the French 
pipeline. SNLS lightcu rves, originally colle cted in qriz 
(Reg nault et al.l [200 9b). are fc-corrected (|Hsiao et al.l 
2007), and standardized by applying a stretch factor, 
to broaden or narrow the rest-fr ame timescale of the 
lightcurve (jPerlmutter et al.lll997| ). in order to generate 
rest-frame B and V band lightcurves3- We define the 
variable r as 

s(l + z) 

where i max is the date of maximum flux (in rest-frame 
B filter band), z is the SN redshift and s the stretch; 
r represents the rest-frame, stretch-corrected, time to 
peak B luminosity. We processed the l ightcurves using 
the SiFTO method (jConlev et al.ll2008D and we used a 
single template to fit the data and stretch correct the 
lightcurves. In processing the SN data we assumed a 
rise time r r = 17.4 days, the time elapsed between ex- 
plosion and maximum B luminosity, according to what 
IGonzalez-Gaitan et al.l (|2011[) finds in a similar (but 
larger) SNLS datas et, and it is al s o cons istent with the 
rise time derived in lHavden et al.l ()2010bl ) from SDSS-II 
SNe la (r r = 17.38 ± 0.17 days). 

Our primary analysis is focused on spectroscopically 
confirmed SNe la. The first 3 years of SNLS data offer 
over 200 spectroscopically confirmed SNe la lightcurves. 
The original dataset was reduced to 87 SN lightcurves 
by applying quality and redshift cuts described below. 
The final dataset uses only SNe that satisfy the following 
requirements: 

• spectroscopically confirmed Type la SNe at redshift 
z < 0.7 (135 lightcurves) 

1 Note that this is different from what is done in the processing of 
SNIa lightcurves for cosmology, where the stretch co rrection is ap- 
plied to the rest-frame tem plate, to match the data llUonlev et alj 
I2U1H ; [Sullivan et alJlMilTl . 
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• the total reduced \ 2 f° r the SiFTO template fit, 
applied to epochs r > — 10 days, is better than 
3.0 (130 lightcurves) 

• the error in the determination of the peak date is 
Ad max < 0.7 days (117 lightcurves) 

• have at least three data points in rest-frame B and 
three data points in rest-frame V band in the rise 
portion of the lightcurve, — 10 < r < days, to en- 
sure the quality of the pre-peak fit (87 lightcurves). 

The excess due to shocking may be visible up to 10 
days after explosion (K10), or r ~ —8 given our choice 
of 7> = 17.4 days. No data prior to r = — 10 days were 
used to standardize the data and generate our composite 
lightcurves in order to avoid including in the lightcurve 
fitting process data points potentially affected by the ex- 
cess that we are seeking. Different choices of minimum 
day (between -10 and -7) were also tested and they do 
not affect our result. However removing points earlier 
than r = —8 causes, in a few cases, a poor lightcurve 
fit, and thus a larger scatter in the data. Visual inspec- 
tion reveals that none of those SNe la for which the fit 
parameters significantly change if data points between 
r = — 7 and —10 are excluded is actually affected by 
shocking. Thus we conclude that including points at 
—8 < r < — 10 only strengthens the significance of 
our results. 

The templates adopted to process the data (Conley09c 
and Conley09fQ, Conlcv ct al. 2008) assume a parabolic 
be havior in time prior to r = — 10 dayfl as d escribed 
in iGoldhaber et al.l (|200l[ ). iConlev et all (f200l . and in 
the references therein : 

f{t)=a{T- Tr )\ (2) 

where / is the flux as a function of time t, and rise 
time r r . This is consistent with a simple expand- 
ing fire ball modeling of t he exploding ejecta (see for 
example Rie ss et ail ll999T ). with a representing the 
rise "speed," which is what we expect in absence of 
shocking by a companion. We refer the reader to 
IGonzalez-Gaitan et al.l ([201 ID for a detailed study of the 
ri se behavior of SNe la in the SNLS data. 

IHavden et al.l ((201 0b) found that a better fit to the 
SDSS-II data can be achieved using two separate tem- 
plates to fit the rise and fall portion of the lightcurves, 
thus obtaining two stretch values. While using 2 
stretches slightly improves the \ 2 P er degree of freedom 
the individual lightcurve fits, F-tests show that for these 
SNLS data the improvem ent achieved using 2 stretche s 
is not significant (see also IGonza lez-Gait an et aLll2011f ). 
Furthermore, because we use only data points at r > 
— 10, with the 5 day cadence of the SNLS data, we gen- 
erally have only 2-3 points in the rise portion of each 
lightcurve that would be used for fitting. Fitting sep- 
arately the rise and fall portions of the lightcurve then 
exposes us to the risk of misfitting or over-fitting the rise 

2 We see no difference in our results choosing either Conley09c 
or Conley09f , and where not specified we will refer to Conley09f 
throughout the paper. 

3 In fact, the SiFTO method allows as well for a cubic fit to 
the early rise portion of the lightcurves, but this was found not to 
improve the lightcurve fit in most cases {Conlcv ct al .1 120081 ). 



portion. We conclude it is best to use a single stretch 
template for the purpose of this analysis. 

A more detailed description of the st andard SNLS 
lightcu r ve processi n g can be fo und in IConlev et al.l 
(|20T1 . IGuv et all d2010h . and IConlev et al.l (120061) . 
For a discussion on the S NLS photometric calibration 
see lRegnault et alj (|2009aft . 

The SNLS lightcurves, normalized to peak flux / pca k = 
1 in each color channel, stretched, and ^-corrected as 
described above, can be combined into a composite 
lightcurve: our composite rest-frame B (V) lightcurve 
contains a total of 1059 (1125) data points between 
r = -20 and r = 40, and 202 (217) in the 10 days 
after explosion that would be affected by the flux ex- 
cess: — 17.4 < t < —7.4 days. The composite B and 
V lightcurves are shown in Figure [TJ B band flux on the 
left hand side and V band flux on the right hand sid^S 
The data points potentially affected by the shocking ex- 
cess are plotted in red, and included within vertical lines. 

3. MODELS 

After an explosion, the SN ejecta expands and collides 
with the companion star, if one exists. K10 showed that 
this impact shocks the SN ejecta creating a hole in the 
expanding material. Radiation can now escape from the 
otherwise optically thick ejecta shell. An early X-ray 
emis sion, analogous to the X-ray flash in core -collapse 
SN (|Soderberg et al.ll200l IModiaz et al.l 120091) . should 
last minutes to hours, with little chance to be observed. 
In UV and optical bands the flux excess lasts longer: 
the gas begins expanding to refill the hole carved by the 
companion. Radiation continues to diffuse out of this 
hot, shocked region, and it is observable until the 56 Ni 
luminosity begins to dominate the lightcurve. The time 
scale for this process is roughly 5-10 days. The effect 
is more prominent in bluer bands and can span over an 
order of magnitude in flux in UV, generating an early 
peak even brighter than the peak luminosity in absence 
of shocking, while it is substantially dimmed in V band. 

The size of the hole and of the shocked gas region, 
and thus the amount of radiation escaping, depend es- 
sentially on the solid angle subtended by the companion 
carving the hole, and on velocity of the ejecta at the time 
of impact. The models assume that the companion star 
is in Roche lobe overflow and that the exploding WD 
has reached the Chandrasekhar limit. Under such as- 
sumptions these are determined by the distance to the 
companion, and the geometry of the system (semi-major 
axis of the binary orbit and size of the Roche lobe) is 
set by the nature of the companion. This excess radia- 
tion is then a powerful tool to identify the SN progenitor 
system. 

K10 considered 3 types of companions: a 2 M Q MS 
star, at a distance a — 5x 10 11 cm from the WD, a 6 M Q 
MS companion at a = 2 x 10 12 cm, and a 1 M RG com- 
panion at a = 2 x 10 13 cm. The effect is most prominent 
for observers viewing directly into, or at small angles 
from, the companion (i.e. into the hole), or roughly 10% 
of the time (Figure [2|). However, particularly in the case 
of RG companions, where the excess is maximum, some 

4 For a discussion on the distinct i on bet ween photon-based, and 
energy-based flux see INugent et all pOOl - ). Throughout the paper 
we refer to flux as photon-based flux. 
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Fig. 1. — Composite rest- frame B (left) and V (right) SNLS lightcurves. Vertical lines delimit the region that could be affected by 
shocking, and the points within are plotted in red. Those are the data points potentially susceptible to the excess induced by the presence 
of a companion shocking the ejecta. 



flux excess is observable even at large angles. We expect 
the distribution of viewing angles to be uniform in na- 
ture, and in our set of lightcurves. The presence of a flux 
excess may increase the detectability of a SN, thus the 
orientation of the WD-companion system may constitute 
a selection bias when operating close to a survey detec- 
tion limit. The cut at z < 0.7 assures that we keep safely 
away from the detection limit of SNLS. Furthermore the 
maximum excess in B band is significantly smaller (a fac- 
tor > 2) than the peak luminosity, and since SNLS is a 
rolling survey, covering each field every five days, each 
lightcurve should, in principle, contain > 3 data points 
within 10 days of peak that would be brighter than the 
maximum shocking flux. The shock-induced flux excess, 
therefore, does not significantly increase the detectability 
of a SN in our sample, and we can assume that the dis- 
tribution of angles in our data is unbiased. Later we will 
expand our analysis to a photometrically selected sample, 
to assess weather an early flux excess might have caused 
the SNe la to be misclassified as a non-SN la, and thus 
not followed spectroscopically (Section [5]) . 

The models are generated from 2-dimensional Monte 
Carlo radiation transport simulations that are subject 
to random sampling errors. Such errors are purely sta- 
tistical, and don't take into account any of the possible 
systematic errors or uncertainties in the model calcula- 
tions. The statistical errors are accounted for throughout 
our analysis. 

The excess generated by shocking is shown in Figure 
averaged over viewing angles, for all three progenitor sce- 
narios in both B and V band. It is evident that, while the 
RG progenitors generate a significant excess, and a very 
distinct effect in both B and V band, the time behavior 
for MS stars is only marginally changed in the presence 
of shocking, especially after averaging over viewing an- 
gles. Such a small deviation from the parabolic early rise 
behavior would hardly be detectable in the presence of 
the typical noise of SNLS data. Therefore we restrict 
ourselves to the RG scenario and only try to constrain 
the RG contribution to SN progenitor systems. We also 
expect that the explosion in a DD scenario would show 
even smaller, or no, deviations from a parabolic behav- 
ior. We thus tentatively associate the DD scenario to the 
standard template. Note, howe ver, that it is possi ble, as 
mentioned in K10 and shown in lFrver et al.1 (|2010t h that 



in a WD merger gas would be blown out to large radii 
(~ 10 13 cm), producing a shock signature, with a UV 
excess possibly propagating t hrough visibl e wave lengths. 
However, the simulations in iFrver et al.l (|2010l ) gener- 
ally produce lightcurves rather dissimilar from standard 
SNe la, with broader visible band lightcurve, unlikely to 
match SNe la in our sample. 

Where needed, we will assume that any explosion not 
generated by a RG-WD binary pair has equal probability 
of arising from any of the three remaining scenarios. 

3.1. Rendering of the models 

The K10 simulations generate full spectra of the SN 
explosion, including the effects of shocking, at time in- 
tervals of 0.1 days for the first 10 days after explosion. 
The spectra are integrated on a day time scale and fil- 
tered through the same V and B filters into which the 
SNLS data have been converted. This is done for every 
angular separation between the WD and the compan- 
ion, in 40 equal intervals of observing probability, for the 
three progenitor scenarios considered: RG, 6 M Q and 
2 Mq MS sub-giants. 

The K10 spectra are designed to reproduce the excess 
generated by shocking on top of a nominal template. 
The input template in the models is irrelevant to the 
shocking physics. We use the spectra for the smallest 
companion scenario {M = 2 Mq) at the largest angu- 
lar separation (~ 180°), where we expect the effects of 
shocking to be entirely negligible, as a neutral template: 
the template in absence of shocking. To better repro- 
duce what we actually expect the result of shocking to 
look like in a SNLS lightcurve we subtract the neutral 
template from the lightcurve templates generated as de- 
scribed above. The new lightcurves are shown in Fig- 
ure [3j averaged over all angles, and they describe the 
excess due to shocking. This excess can be added to the 
parabolic portion of the templates to reproduce what 
we would expect to see in our data. The reader is re- 
minded that this portion of the lightcurve is not used to 
standardize the SNLS data and generate the composite 
lightcurves. 

We now have template lightcurves for the first 10 days 
of a SN la explosion in the SD scenario, with different 
companion stars and at different observing angles, which 
can be compared to the SNLS observations. 
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Fig. 2. — The K10 model for a WD accreting from a RG companion is shown. At the center is a schematic representation of the SD 
explosion scenario: in the expanding ejecta, gray, the impact with the companion star (black circle) has created a hole, here simplistically 
represented by a cone. To the left and right, according to the corresponding view point, are the rise lightcurves for, respectively, an 
observer looking in opposite direction from the companion (no excess), and looking into the companion and the hole created by the impact 
(maximum excess), for the case of a WD-RG progenitor system. The scatter in the model is simply due to statistical noise in the numerical 
simulations (see Section [3}. The solid line is the Conley09f template. 
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Fig. 3. — Models of excess emission over a nominal - parabolic - SN lightcurve template, in units of peak luminosity, signature of shocking 
by a companion star, for the cases considered in K10: a RG companion (red filled circles), a 6 Mq (blue empty circles) and a 2 Mq (green 
crosses) MS companion, all in Roche lobe overflow (separation from the core of the explosion a = 5 X 10 11 , 2 X 10 12 , and 2 X 10 13 cm 
respecti vely ). The effect is shown averaged over all observing angles. SN la spectra are generated from KlO's simulations as described in 
Section |3. II and filtered through standard B and V filters to generate the theoretical lightcurves. The error bars represent the scatter - 
standard deviation - in the models. The left plot shows the effect in B and the right hand plot in V band. 



4. TESTS 
4.1. Template goodness of fit 

We begin by noticing that the fit of the composite SN 
to the SN template (here Conley09c and Conley09f were 
used, with no significant differences) is worst in the re- 
gion of interest for the shocking effect, r < — 7 days to 
peak, in both B and V bands. 

Figure U shows the median, binned by day, of the com- 
posite lightcurves (top plots). The error bars represent 
the scatter - standard deviation - of the individual mea- 
surements (the standard deviation is measured as con- 
ventionally done with respect to the mean and we ignore 
the error of each measurement). The bottom plots show 



the deviation from the template as the difference between 
the data, /, and in the template, T, at time t, over the 
error in the data o~(t), averaged over each day: 



E 

t in day 



/(*) ~ T(t) 
a{t) ' 



(3) 



as an estimator of the goodness of fit of the template 
to the composite lightcurve. The propagation of errors 
along the time dimension is also ignored. 

The most significant deviation from the template hap- 
pens in both B and V band roughly prior to r = — 10 
days, with a clear excess in B band at r < —12 or in 
the first few days after explosion. It is intriguing that the 
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Fig. 4. — Top: median of the composite lightcurve shown in Figure [T] in rest-frame B (left) and V (right). The error bars represent the 
scatter in the measurements in the composite lightcurve. Bottom: a measure of the deviation of the median plotted in the top panels from 
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the nominal SN templates, averaged day by day: 
t =-11 days in B and before r =-7 days in V 

deviation is more prominent in B than in V band, con- 
sistent with the chromatically biased effect that shocking 
by a companion would produce. Note, however, that this 
is the portion of the lightcurve that is not fit to the tem- 
plate, and it is thus not surprising to see a larger scatter 
here. 

4.2. Simulations 

Having found a deviation from our fiducial SN la tem- 
plate in the early days after explosion, we test if this can 
be attributed to shocking by a companion. 

With the K10 templates in hand we can create syn- 
thetic SN la rise lightcurves that incorporate the effect 
of shocking for the different progenitor scenarios, as seen 
from different viewing angle. We start off with the stan- 
dard parabolic rise templates (Conley09f ). In each band 
we add the excess described by the models (Figure [3} to 
our standard template. For each epoch corresponding to 
the SNLS data, we draw a data point from the new tem- 
plate thus obtained. To choose the viewing angle from 
which this data point should come we draw angles with 
equal probability between and 180°. 

Families of synthetic lightcurves are generated using 
increasing contributions of data points from the RG 
progenitor scenario, and drawing the remaining points 
equally likely from a parabolic template (DD scenario), 
the 2 Mq and the 6 M Q MS progenitor scenarios. We 
generate families with 0% RG contribution, to 100% RG 
contribution, in steps of 10%. A finer grid was also 
tested, but a resolution of 10% in the RG contribution is 
adequate to represent the progenitor populations given 
our errors. 

In our test we have to account for the errors in both 
the templates and the data. Thus at a given epoch, in 
generating the simulated data points, we draw each flux 
value from a Gaussian distribution around the template 



The deviation of the composite lightcurves is most significant before 

value at the corresponding epoch, where the width of the 
Gaussian is the sum in quadrature of the errors in the 
SNLS data and in the model. 

For each RGf rac between RGf rac = 0% and 100% we 
create 100 sets of simulated observations, in steps of 10%, 
thus generating 100 synthetic lightcurves per RGf rac , 
each one the size of the rise portion of the lightcurve: 
202 points in B and 217 in V band. 

4.3. One band K-S test 

We then compare each population of synthetic 
lightcurves to our composite lightcurves. We chose 
the non -parametric 2- sample Kolmogorov-Smirnov (K- 
S) test (|Peacocldll983f ) to do so. This simple statistical 
test measures the maximum distance between the simu- 
lated and true cumulative lightcurves, and we find it is 
a sensible statistic to determine the presence of an ex- 
cess in data with large scatter. The K-S tests are non 
parametric, and thus do not require us to make any as- 
sumptions on the distribution of d ata, and, u nlike for 
example the Pearson's \ 2 statistics (Rice 2001), do not 
require binning the data. The SNLS data points are 
sorted by flux, and added to generate a cumulative flux 
distribution. Similarly, a cumulative flux distribution is 
generated for each synthetic realization. The maximum 
distance between two cumulative distributions is a mea- 
sure of the level at which the null hypothesis that the 
data sets being compared come from the same distri- 
bution, can be rejected (see Figure [5]). In other words: 
comparing the true data with synthetic data generated 
from scenarios with different RG contributions, we test if 
the data comes from a distribution with a given fraction 
of RG progenitors, RGf rac . 

For each RGf rac , we pair each of the 100 sets of sim- 
ulated observations with our true lightcurve, apply a 2- 
sample K-S tests, and average over the 100 K-S num- 
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Fig. 5. — Cumulative, normalized flux distribution for the rise 
portion of the SNLS composite B band lightcurve (solid line) and 
for one of the synthetic realizations (dotted line), with a contribu- 
tion of RG progenitors RGf rac = 30% (see Section l4.2H . The data 
points in the lightcurves are ordered by flux and the x-axis is simply 
the rank of each sorted data point. The maximum vertical devia- 
tion between the true and synthetic cumulative flux distributions 
is marked by an arrow. From this value we asses the probability 
that the true data and this synthetic realization came from the 
same original distribution. 



bers. This allows us to obtain the confidence level for 
the rejection of the null hypothesis - that the two sets 
being compared come from the same distribution - and 
its statistical errors, accounting for all sources of noise: 
the uncertainty in the models, in the data, and for the 
presumed diversity in SN la progenitors. 

Using only the B data, and within la error bars, we 
reject at roughly 90% confidence level (c.l.) the hy- 
pothesis that a progenitor population with RGf lac > 
30% has generated our data. A contribution of more 
than 40% is ruled out at the > 95% c.l. These re- 
sults are large l y consi stent with the analysis presented in 
lHayden et al.l (l2010af). and w ith the upper limits placed 
bv (|Gilfanov fc Bogdanll20i0l) in SNIa in elliptical galax- 
ies. Accounting for all angles, with only ~ 200 data 
points in the 10 days that would be affected by the shock- 
ing effect, we would expect < 40 data points to be af- 
fected significantly by the effect we are probing, even 
if all companions were RGs. This is therefore a small 
number statistics problem, in the presence of noise in 
both data and models, and it is not surprising that we 
have a limited ability to place strict limits to the con- 
tribution of RG progenitors this way. However we will 
obtain stronger limits by including considerations on the 
color bias in the shocking excess in the next section (Sec- 
tion H3). 

Noticeably, despite the noise, the probability that the 
true and synthetic distribution of data points come from 
the same progenitor distribution clearly decreases mono- 
tonically as we increase the contribution of RG systems, 
particularly in the B band (Figure |6|). This strongly 
suggests a minimal contribution of RGs to SN la pro- 
genitors. Similarly, an almost monotonic decrease in 
probability (increase in c.l.) is evident in the V band, 
though less pronounced. This is expected, on account 
of a smaller signature of shocking in redder bands (see 
Figure EJ. 
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Fig. 6. — Results from 2-sample K-S tests applied to the early 
rise portion of the composite lightcurves from the SNLS. We plot 
the confidence level (c.l.) at which we can reject the null hypothesis 
that the true data comes from the same distribution of progenitor 
that generated the simulated data, as a function of RGf rac . The 
90% c.l. is marked with a solid black line. The true and synthetic 
population are increasingly less likely to come from the same dis- 
tribution as the RG contribution to the mixture of synthetic pro- 
genitors is increased: RGs are unlikely to be common progenitors 
of SNe la. The effect is more pronounced - and strictly monotonic 
— in B band (blue) than in V band (red) because of the color bias 
in the shocking excess. A RGf rac < 0.3 is allowed within the ltr 
error bar while larger contributions are ruled out at 90% c.l. or 
greater 

4.4. K-S chromatic test 

In this section we investigate the chromatic bias in the 
shocking footprint. In absence of shocking, the expected 
time behavior of the SN explosion is a parabola, simi- 
lar in V and B band. Thus, in the rise portion of the 
lightcurves we would expect points drawn from a set of 
SNe to come, statistically speaking, from the same dis- 
tribution in V and B band. However, in the K10 sim- 
ulations (see Section [3]) the V and B time behavior dif- 
fer dramatically in the rise portion of the lightcurve in 
the presence of RG progenitors. We again perform a 
2-sample K-S test. This time we want to asses the sim- 
ilarity of the B and V populations of early-rise data, so 
for the SNLS data and we compare the B and V channel 
with a K-S test, and we do the same for each synthetic 
population. 

We find that the hypothesis that rest-frame B and 
V populations of data points from the composite true 
SNLS lightcurves, day 0-10 after explosion, come statis- 
tically from the same distribution can only be rejected 
to < 5% c.l., or equivalently that the hypothesis that 
the two channels come from the same distribution has a 
p — value ~ 0.95. 

We compare the synthetic B and V lightcurves, and 
find, as expected, that the K-S number increases with 
the increasing RG contribution: the probability that the 
B and V synthetic data come from the same distribu- 
tion decreases as more RG progenitors are used in the 
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Fig. 7. — Results from 2-sample K-S tests performed to compare 
the early rise portion of the V and B band lightcurves. The solid 
line represents the p — value of the hypothesis that the true V and 
B lightcurves come from the same distribution: p — value ~ 0.95; 
below the line the B and V would diverge more than the SNLS 
rest-frame B and V do. The points represent the p — value's for 
our synthetic B and V data, as a function of the RG contribution 
in the progenitor mixture. That is: p — value of the hypothesis 
that, for each RG contribution, the B and V band synthetic data 
came from the same distribution (where the p — value is 1 minus 
the rejection level of the hypothesis). The la error bars plotted are 
obtained generating 100 synthetic populations and accounting for 
the errors in the SNLS data and in the K10 simulations. Only the 
synthetic population that contains no RG progenitors is consistent 
at the lcr level with the data. A RGf rac > 10% is ruled out at the 
2<r, and a RGf rac > 20% at the 3cr level. 

progenitor mix. 

The results of this test are plotted in Figure [71 Within 
la error bars, only the population with no contribution 
from RG progenitors is consistent with the SNLS data. 
We rule out a contribution RGf rac > 10% at ~ 2a, and 
greater than 20% at > 3a level. 

4.5. Color distributions 

In standardizing our lightcurves we have chosen one 
stretch value for each lightcurve to be applied to all rest- 
frame bands. We ask the question: could there be a 
correlation between V and B that would interfere with 
the result of our K-S chromatic test (see Section I4.4[) . 
Suppose some data points in the region r < — 10 days, 
which is not fit to the template, have more (or less) flux 
than the template, so that if we included those points 
when fitting the lightcurve to the template we would have 
generated different fit parameters, and a different value 
for the stretch or day of maximum; in this case both the 
B and V lightcurves would show flux in excess (deficit) 
of the template in the shock region. We might then see 
a correlation in our B and V composite lightcurves. Any 
such correlation would systematically lower the K-S num- 
ber (higher the p — value for the null hypothesis) found 
for the true data. Meanwhile our simulated lightcurves 
are generated independently in B and V bands. Since in 
our chromatic K-S test even the simulations which use no 
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Fig. 8. — fB— v distribution for true (top) and simulated data 
(four bottom plots). The synthetic populations are generated 
from progenitor mixtures with, from the second to the last plot, 
R.G frac =0, RG frac =33%, RG frac =66%, and RG frac =99%. The 
mean for each distribution is plotted as a vertical arrow, and the 
standard deviation of each population as a double-ended arrow 
above the distribution. A dotted-dashed line shows the location of 
fe-v =0. B and V are normalized at peak, thus no color excess 
is expected in absence of shocking. Each synthetic population is a 
factor 100 larger than the true population, to minimize the Poisson 
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data at the la level, but show a lower p — value, indi- 
cating a greater discrepancy between B and V band, we 
further explore this possibility, seeking an independent 
confirmation of what we see in the K-S tests beyond this 
possible systematic correlation. 

We test the color of our true and simulated lightcurves 
by taking fe-v to be the difference of the B and V flux 
after normalizing each channel at peak. Under the as- 
sumption that in absence of shocking the rise portion of 
the lightcurve would follow a parabolic behavior identical 
in both bands, diverging only at r > —9 as modeled in 
the Conley09f template, the distribution of fe-v values 
should be consistent with for our composite lightcurves, 
while the effect of shocking would produce a distribution 
of fe-v with a positive mean, and a large standard de- 
viation. 

For every flux point in each normalized, standardized 
rest-frame B lightcurve, we subtract the flux of the clos- 
est rest-frame V data point, within A r < 0.2 days. Sim- 
ilarly, we generated synthetic colors for different RG frac- 
tions, by creating pairs of V and B synthetic lightcurves, 
accounting as usual for the typical error bars in the data 
and in the models. We derive the distribution of colors 
for both true and simulated data. 

The fs-v distributions are plotted in Figure[5] The top 
panel shows the fs-v distribution for true data. There 
is no blue excess in flux in the true color: in fact the 
distribution has a mean of ^ ~ — 8 x 10~ 4 , a median ~ 
0.002, and a standard deviation a ~ 0.054: statistically 
consistent with a random distribution around 0. 

The distributions generated from simulated 
lightcurves are shown below the distribution 
for true data in Figure HJ f° r RG contributions 
RG frac = 0%, 33%, 66%, and 99%, plotted from 
the top to the bottom. Each distribution is generated 
from a factor of 100 more points than the true color 
distribution and is thus minimally noisy. The mean of 
the distribution increases as we increase RGf rac and the 
distributions get increasingly asymmetric, weighted to- 
ward positive values of fe-v (bluer color). The synthetic 
distributions generated with no RGs (RGf rac =0%) has 
moments that are extremely similar to those of the true 
color distribution: (i ~ 2x 10~ 3 , median ~ 0.001 and 
a ~ 0.077. Once again, this shows that the distribution 
of colors in the SNLS data is compatible with minimal 
- or no - contribution of RG to SN la progenitors, 
confirming the results obtained from the K-S tests. 

5. PHOTOMETRICALLY SELECTED SNe la 

The excess due to shocking of the SN ejecta affects the 
early time domain photometric and spectral behavior of 
the SN la explosions. Since in surveys such as SNLS and 
SDSS SN la are identified by their early lightcurves, and 
thus an explosion is followed up spectroscopically only if 
it is thought to be a SN explosion, an interesting ques- 
tion is whether this early effect might have lead to the 
rejection of phenomena that indeed were SN la, but de- 
viated from the expected early be havior on account of 
shocking. In Hav den et al.l (|2010al ). a subset of uncon- 
firmed SN la is visually inspected and no such effect is 
found. We investigate 905 SNLS lightcurves with some 
rcdshift information, either spectroscopic or photomet- 
ric. We exclude likely or known AGN, variable stars, 
and core-collapse (CC) SNe. In order to avoid contami- 



nation from unidentified SNe II, lb, or Ic, we also apply 
cuts in stretch and color space. In particular, CC SNe 
show a different average color than SNe la, and color 
constraints eliminate them from the sample. A detailed 
discussion of photo metric selection of SN e la in the SNLS 
data can be found in lBazin et al.1 ()2011h . We thus believe 
our new dataset has minimal contamination from non 
SN la events. Our new dataset contains 336 lightcurves 
before our cuts are applied (see Section [2]), and 110 af- 
ter. Our new composite lightcurves contain 251 points 
in rest-frame B and 270 in rest-frame V in the region of 
interest: r —-17 A to -7.4 days to peak (Figured]). 

We repeat the K-S tests applied earlier to the extended 
SN la set and find that the statistics confirm the upper 
limits set to the contribution of RG binary systems to 
SN la explosions (Figures Hfl] and [TT])- The K-S test 
of the composite lightcurve in each B and V with the 
respective synthetic lightcurves is entirely consistent with 
the test for the spectroscopically confirmed SN la subset, 
and consistent with minimal or no contribution of RG to 
the SN la progenitors. 

6. U BAND DATA 

As described in Section [31 the excess due to shocking 
is more prominent at bluer wavelengths. In U band, the 
models by K10 predict an excess over a nominal template 
more pronounced by roughly 20% over the B band for the 
RG case, averaged over all angles and over the first 10 
days after explosion. The prediction for the angle aver- 
aged excess in U band is shown in Figure [121 We thus 
extend our analysis to the rest-frame U band data, to see 
if a stronger constraint can be placed to the contribution 
of RG to the SN la progenitor population. Three years 
of SNLS spectroscopically confirmed lightcurves are pro- 
cessed as described in Section [21 in order to generate 
rest-frame U band lightcurves. The lightcurves are then 
selected if they pass similar cuts to those described ear- 
lier: 

• spectroscopically confirmed Type la SNe at redshift 
z < 0.7 (135 lightcurves) 

• the total reduced % 2 for the SiFTO template fit, 
applied to epochs t > — 10 days, is better than 
3.0 (130 lightcurves) 

• the error in the determination of the peak date is 
Ac? max < 0.7 days (117 lightcurves) 

• have at least three data points in rest- frame B, 
three data points in rest-frame V, and three in rest- 
frame U band in the rise portion of the lightcurve, 
— 10 < t < days, to ensure the quality of the 
pre-peak fit (57 lightcurves). 

The rest-frame U band composite lightcurve thus gen- 
erated is plotted in Figure Q21 and it contains 662 points 
between days -20 and 40 from explosion, and 123 in the 
region of interest: — 17.4 < t < — 7.4. Applying the lat- 
ter cut, which is more restrictive than the corresponding 
cut in our primary analysis, the new composite B and 
V lightcurves contain, respectively, 152 and 161 data 
points in the region —17.4 < r < —7.4. It is imme- 
diately evident that the U band composite lightcurve is 
significantly noisier than the B and V composites, and it 
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Fig. 9. — Composite B (left) and V (right) SNLS lightcurves containing both confirmed and unconfirmed SNe la. Points plotted in red, 
and enclosed by vertical lines, are those which may be affected by the excess induced by the presence of a companion. 
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Fig. 10. — Results from 2-dimensional K-S applied to the early 
rise portion of the composite lightcurves. The Figure reproduces 
Figure [6] for the extended SN la set, not limited to spectroscopi- 
cally confirmed events. 



Fig. 11. — Results from 2-dimensional K-S applied to the early 
rise portion of the color composite lightcurves. This Figure repro- 
duces Figure for the extended SN la set, not limited to spcctro- 
scopically confirmed events. 



contains 35% fewer lightcurves, and roughly 40% fewer 
relevant data points. 

We generate synthetic U band lightcurves as described 
in Section |4~21 and we reproduce the 2-sample, monochro- 
matic K-S test we described in Section 14731 for the U band 
data. We find that the true data distribution once 
again grows dissimilar from the simulated distribution 
as more RG progenitors are included in the simulation 
(Figure fl4"|) . Using the U data a progenitor fraction 
RGf rac >= 30% is ruled out at the 3cr level. The re- 
sults of the single band K-S test for the B and V for the 
subsamplc of lightcurves that pass the new set of cuts 
arc also plotted, and they are entirely consistent, though 
with larger errors on account of the smaller dataset size, 
with the results obtained in Section 14.31 

Note that the U band data, even in absence of RG in 
the simulated data, appears from a K-S test to be dif- 
ferent at the 2<T level from the parabolic-rise model (i.e. 



the 2cr limit of the c.l. of rejection of the null hypothesis 
that synthetic and true data come from the same distri- 
bution is below c.l. = for all synthesized populations). 
The simulated lightcurves are generated as described in 
Section 14.21 and an adiabatic (parabolic) expansion is 
postulated up to 6 days after explosion. However, effects 
of line opacity and dispersion in the spectr a at wave- 
length s bluer of A = 400 nm, as described in lEllis et aTl 
(2008), affect the U lightcurve, and may modify it from 
the our simple model prediction. This can explain the 
relative low correlation between our true and simulated 
data, which is revealed by the K-S test. For this reason 
we are reluctant to apply the chromatic K-S test describe 
in 14.51 to the U band data, as our assumption that the 
U and B, or U and V lightcurves would have identical 
early rise behavior might not hold here. As a matter 
of visualization though, we reproduce Figure [8] for the 
fu-v an d fu-B data (see Section |4~5|) . As in all other 
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Fig. 12. — Excess over a template caused by shocking in U rest- 
frame band, according the the theoretical models presented in K10: 
the figure reproduces Figure [3] for the U band. The angle averaged 
excess, in units of peak luminosity, is shown for a RG with red 
filled circles, for a 6 Mq with blue empty circles and for a 2 Mq 
with green crosses. 
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days to peak 

Fig. 13. — Composite U lightcurve from three years of SNLS 
data. Symbols are as described in Figure [T] 

cases considered, the distribution of simulated color gets 
increasingly different from the distribution of true data 
as we increase the contribution of RG in the progenitor 
mix for our simulations. 

We limit ourselves to point out that the U band data 
confirms the constraints that we set with B and V data. 

7. CONCLUSIONS 

We analyzed 3 years of spectroscopically confirmed 
SNe la from the SNLS survey looking for an early rise 
flux excess that could be attributed to shocking by a 
companion. We created composite lightcurve standard- 
izing the data with the SiFTO method, excluding from 
the fit the region that might be affected by the shocking 
phenomenon. 

We found a worsening of the fit of the data to the 
standard templates in the first few days after explosion 
(Section 14. ip . but we found no evidence that is due to 
anything but the fact that the data prior to r = — 10 
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Fig. 14. — Results from 2-dimcnsional K-S applied to the early 
rise portion of the composite U (black crosses), B (blue full circles), 
and V (red empty circles) lightcurves; la error bars are shown. 
This Figure reproduces Figure \E\ in the three color bands, and for 
the subset of lightcurves selected by the cuts described in Section|6] 
In U band, a RG fraction of 30% and its 3a error bars lie below the 
solid line, which indicates the 0.95 c.l of rejection of the hypothesis 
that the simulated and true data come from the same distribution. 

are not used in the template fit. 

We used the spectra generated by the K10 (Kasen 
2010) simulations to model the expected SN la time do- 
main behavior in the SD scenario, thus we can account 
for sources of noise in the models, as well as in the data. 
We found no evidence of flux excess in our data, and 
conclude that, based on the K10 models, the contribu- 
tion from RG progenitors is less than 10% in the SNLS 
3-year sample. We thus set a ~ 2a upper limit of 10% 
to the contribution of RG-WD binary systems to the 
SN la progenitors, and a 3<r upper limit of 20%. With 
roughly 100 lightcurves in our sample, with a contribu- 
tion of ~ 10% lightcurves from RG progenitors, ~ 3 data 
points could be affected by shocking. We cannot exclude 
such a small contribution from RG binary systems in the 
presence of noise from both the models and the data. 
Our results are robust when tested in a photometrically 
selected sample of lightcurves, as well as using U band 
data. 

Our conclusion agre es with the results derived 
in lHavden et all (|2010af ) from the SDSS-II SN la sam- 
ple. Our analysis differs, other than in the SN sam- 
ple u sed, in the treatment of the shocking signature: 
while iBavden et al.l (|2010a) models the shocking as a 
Gaussian excess we used the K10 simulation directly 
to characterize effects of shocking, thus including the 
uncertainties in the models. Furthermore, our analy- 
sis exploited the color bias in the shocking excess to set 
stronger constraints on the presence of shocking and are 
able to quantify the maximum allowed contribution of 
RGs to the SN la progenitors. 

Although Bayesian tests (Expectation Minimization 
and Gibbs sampling) were applied to our data, the pres- 
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Fig. 15. — fu— v (left column) and fu— B (right column) distributions for true (top) and simulated data (four bottom plots) with different 
contribution of RG progenitors. The mean and standard deviation of each population are indicated by arrows. This Figure reproduces 
Figure 151 for the frj— V an d fu— B colors. 
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ence of noise, and the relatively high dimensionality of 
the problem, with four possible progenitor scenarios, 
does not allow us to firmly asses what contribution of 
RGs best reproduces the residuals we see in the data 
with respect to the parabolic templates. Our data is en- 
tirely consistent with no RG progenitors. 

Accor ding to population synthesis studies, such as for 
example iRuiter et al.l (|2009t ). the SD scenario is expected 
to produce SN la mainly from evolved companions, i.e. 
they favor the RG-WD channel over the MS- WD chan- 
nel, at l east under the Ro che lobe overflow require- 
ment. In IRuiter et al.l (|2009f ) . where reaching the Chan- 
drasekhar limit is required, as it is assumed in this paper 
and in the K10 simulations, MS- WD binaries are respon- 
sible only for 5% to 10% of the SN la production, while 
the majority of SNe la come from a system with an 
evolved donor: a sub-giant, or giant. Limiting the con- 
tribution of RG-WD SN la progenitors from an observa- 
tional point of view may then have a significant impact 
on the conclusions derived from population synthesis sti- 
dies on delay time distributions and SN la progenitors. 

The PanSTARRS Medium Deep Survey (PS1, 
iPastorello et al.ll2010f) and the Palomar Transient Fac- 
torv (PTF. lLaw et al.ll2009l ) have begun providing well- 
sampled rise lightcurves, where the excess due to shock- 
ing by a RG progenitor, should this be a valuable chan- 
nel to produce SN la, could soon be observed. Since 
the effect is predicted to be chromatically biased (see 
Section [3|), PS1 is particularly suitable, offering data in 
SDSS g and r bands, thus allo wing a color comp ari- 
son. Early UV follow up surveys ()Cooke et al.ll2011l) are 
also a promising way to spot WD-RG progenitor sys- 
tems, since the progenitor excess is extremely prominent 
in UV bands, provided that the follow up can be trig- 
gered early enough after explosion. SNLS continued col- 
lecting SN la time series through 2006 and as all SNLS 
data become available more stringent limits may be set. 
Note that in absence of any excess the population of pro- 
genitors could be pinpointed to small (M < 6 M Q ) 
progenitor companions. However, in the presence of de- 
tections of small excess signals there would be a degen- 
eracy between RG progenitors observed at some angu- 
lar offset from the line of sight to the hole generated by 
the companion, and more massive MS companions, and 
a large sample is indeed necessary to disentangle these 
two scenarios. Possibly, only a su rvey as large as LSST 
(jLSST Science Collaboratidl2009| ) would offer the oppor- 
tunity to asses the frequency of progenitor companion 
types in SNe la. 

We also point out that the constraints derived here rely 
on the theoretical models described in K10. The shock- 
ing signatures predicted in K10 assume the companion is 
in Roche lobe overflow, with the separation distance, a, 
only a few times the stellar radius, R. While this is ex- 
pected in a typical accretion scenario, if a >> R the solid 
angle subtended by the companion would be smaller , and 
so would be the effect of shocking. Uusthaml (|2011l) . for 
example, argues that the donor star in the SD scenario 
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might shrink rapidly before explosion, having exhausted 
its envelope; the companion star would then be many 
times smaller then its Roche lobe, reducing the shock- 
ing signature, and also explaining the lack of hydrogen 
in spectra of SNe la. We look forward to more detailed 
theoretical work, which may relax the Roche lobe over- 
flow assumption, integrate three dimensional explosion 
models, and takes into account possible absorption mech- 
anisms within the systems, and the effects of the orbital 
motion, to better characterize the shocking behavior and 
its diversity. 
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